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CLAIMS 

What is claimed is: 

1. A natural gradient Blind Multi User Detection (BMUD) network 
system that adaptively estimates a set of matrices to counter a linear convolutive 

5 environment model r n , the system comprising: 

an input receptive of at least one of the linear convolutive 
environment model r n or a whitened version r n M ' of the linear convolutive 

environment model r n ; 

parametric matrices W 0 and W k (k=1,2,... K) adaptable to estimate 
10 independent user symbols y n a\ an n th instant based on at least one of the 
linear convolutive environment model r n or the whitened version r„ w of the 
linear convolutive environment model r n ; and 

a decision stage interpreting y n and estimating corresponding user 
symbol estimates b n also at the n th instant. 

15 

2. The system of claim 1, wherein the system is networked in a 
feedforward configuration. 

3. The system of claim 2, further comprising a recovery stage 
20 adapted to compute y n according to: 

y n =w 0 r:+f i w k c k , 

jfc=l 

where K is an estimate of a number of a previous symbols needed for 
computation of y n , with K being greater than or equal to J:= integer (max 
(Tau_L)) + 1). 

25 

4. The system of claim 2, wherein the parametric matrices W 0 and W k 
have update laws according to: 

AW 0 - (/ - yiyjyf* )w 0 ; and 
15 
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AW k - (/ - <p{y n )yH )w k - q*yjg_ k f , 



where is an element-wise acting score function, I is a K-d identity 
matrix, and k=1,2,...K. 

5. The system of claim 2, wherein W Q is initially chosen to be at least 

one of an identity or a diagonally dominantly matrix, while all other matrices W k 

are initialized to have at least one of random elements with a very small variance 
or as matrices of all zeros. 

6. The system of claim 1, wherein the system is networked in a 
feedback configuration. 

7. The system of claim 6, wherein the recovery stage is adapted to 
compute y n according to: 



8. The system of claim 6, wherein the parametric matrices W Q and W k 
have update laws according to: 



where is an element-wise acting score function, I is a K-d identity 

matrix, and k=1,2,...K, with K being an estimate of a number of previous 
symbols needed for computation of the parametric matrices, K being greater 
than or equal to J:= integer (max (Tau_L)) + 1). 

9. The system of claim 1, wherein the system is networked in a 
feedback configuration without need for any matrix inversion. 



y n =Wo 1 r™-zw k y n - k 



V k=l J 




AW k ~W 0 (p(y n )yX), 
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10. The system of claim 9, wherein the decision stage is adapted to 
compute y n according to: 

K 

k=l 

5 11. The system of claim 9, wherein the parametric matrices W 0 and W k 

have update laws according to: 

AW 0 - (/ - <p(y n )yH )w 0 ; and 
AW k - (/ - (p{y n )yH )w k + <p(y n )yH_ k , 

where is an element-wise acting score function, and I is a K-rf identity 
1 0 matrix. 

12. The system of claim 1, further comprising a whitening filter 
preprocessing received data for dimension reduction to K , which is an actual 
number of principal independent symbol sequences in the received data, and to 

15 remove second order dependence among received data samples and additive 
noise. 

13. The system of claim 12, wherein the whitening filter whitens data 
online using adaptive principle component analysis computational techniques. 

20 

14. The system of claim 13, wherein the whitening filter whitens data 
using an algebraic PCA estimate over a large batch of received data including N 
samples according to: 

R = [r x r 2 ... r N-l r Nj 

25 with a data correlation matrix 

A c =— i— RR T . 
c N-l 

15. The system of claim 14, wherein the filter achieves the whitening 
using a filtering matrix according to: 

W = D' in V T , 

17 
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where D represents a K-dim matrix of principle eigenvalues of the data 
correlation matrix A c , and V represents a KxN matrix of principle eigen vectors 

of the data correlation matrix A c , with K representing a number of users. 

5 16. The system of claim 12, wherein the filter is adapted to calculate 

the whitened version r» of the linear convolutive environment model r n 
according to: 

rn = W(H 0 b n + Hib n _ x +n n ) = H 0 b n + #iV-l . 
and the linear convolutive environment model r n is represented according to: 

where b n and b n ^ are the K-d vectors of current and previous symbols for all 
the K users, H 0 and H x are KxK mixing matrices with the structure 

H 0 =|H 0 ,o H oa ... Ho^H^Lh^ H w ... H liK J such that 

L-l L-l 

Ho,k=Tj€oYi h iZu> Hu=V^"2 A /^i and represent the energy of the 

z=o z=o 

15 current and the previous symbol respectively. 

17. An adaptive detector utilizing knowledge utilized by a RAKE 
receiver, comprising: 

an adaptive weighting matrix introduced into a RAKE structure, 
20 wherein the matrix is adaptively estimated using at least one of Principal 
Component Analysis (PCA) computational techniques and static Blind Source 
Recovery (BSR) computational techniques. 

18. The detector of claim 17, further comprising a closed form RAKE 
25 structure according to: 

£ | S&WH H r for DS-CDMA Systems 

URAKE-ICAI PCA " 1 S H C H^h h 7 f or WCDMA Systems 

where W = diag[AA • • • A ] , and A is the matrix. 
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19. The detector of claim 18, wherein W{or A) is adapted according to 
natural gradient update laws. 

20. The detector of claim 19, wherein the matrix is adaptively 
5 estimated using natural gradient update laws according to: 

A(k + 1) = A(k) + TJ k AA(k) . 

21 . The detector of claim 20, where 

AA(Jk) = {( J "^( fe ))yW ff ) A W for static BSR{orICA) 
\(I-y(k)y(k) H )A(k) forPCA 

10 and <p(-) is a nonlinear score function which depends on an underlying 
distribution structure of involved signals. 

22. The detector of claim 21 , wherein the score function is in the form: 
<Pi 0>/ ) = 1% yi ~ oti (tanh(# Refo }) + tanh(# Im{y t })) . 

15 

23. The detector of claim 17, wherein a channel estimate H is at least 
one of not available or changes dynamically, and the detector is estimated 
without using the channel estimate, such that the structure reduces to Matched 
Filter (MF) BSR/PCA according to: 

J StfWr for DS-CDMA Systems 

^ J iMF-BSR I PCA ~ | S P C Hwr forWCDMA Systems 

24. A natural gradient Blind Multi User Detection (BMUD) method that 
adaptively estimates a set of matrices to counter a linear convolutive 
environment model r n , comprising: 

25 receiving at least one of the outputs of the linear convolutive 

environment model r n or a whitened version r w w of the outputs of the linear 

convolutive environment model r n ; 
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adapting parametric matrices W Q and W k to estimate independent 
user symbols y n at an n th instant based on at least one of the linear convolutive 
environment model r n and the whitened version r* of the linear convolutive 
environment model r n ; and 

interpreting y n and estimating corresponding user symbol 
estimates b n also at the n th instant. 

25. The method of claim 24, further comprising employing a 
feedforward network configuration. 

26. The method of claim 25, further comprising computing y n 
according to: 

k=l 

27. The method of claim 25, further comprising updating the parametric 
matrices W 0 and W k via update laws according to: 

AW 0 oc(/^^(3; n )3;^)wo;and 

AW k - (/ - piyJyV )w k - <P(y n ){rZk f - 
where is an element-wise acting score function, and I isa K-rf identity 
matrix. 

28. The method of claim 25, further comprising: 

initializing W 0 to be at least one of an identity or a diagonally 
dominant matrix; and 

initializing all other matrices W k to have at least one of random 
elements with a very small variance or as matrices of all zeros. 
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29. The method of claim 24, further comprising employing a feedback 
network configuration. 

30. The method of claim 29, further comprising computing y n 
according to: 



31. The method of claim 29, updating the parametric matrices W 0 and 
W k via update laws according to: 



where q>{) is an element-wise acting score function, and I is a K-d identity 
matrix. 

32. The method of claim 24, further comprising employing a feedback 
network configuration without need for any matrix inversion. 

33. The method of claim 32, further comprising computing y n 
according to: 



34. The method of claim 32, further comprising updating the parametric 
matrices W 0 and W k via update laws according to: 



V k=l J 




k=l 




&W k ~[l-<p(y n )y%)w k +<p(y n )y 



'n-k ' 
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where is an element-wise acting score function, and I is a K-d identity 
matrix. 

35. The method of claim 24, further comprising preprocessing received 
5 data for dimension reduction to K, which is an actual number of principal 

independent symbol sequences in the received data, and to remove second 
order dependence among received data samples and additive noise. 

36. The method of claim 35, further comprising whitening data online 
10 using adaptive principle component analysis computational techniques. 

37. The method of claim 36, further comprising whitening data using an 
algebraic PCA estimate over a large batch of received data including N samples 
according to: 

15 R = [r l r 2 r N-l r Nj 

with a data correlation matrix 

1 T 
c N-l 

38. The method of claim 36, further comprising employing a filtering 
20 matrix according to: 

W^D~ V2 V T 

where D represents a /C-dim matrix of principle eigenvalues of the data 
correlation matrix A c , and V represents a KxM matrix of principle eigen vectors 

of the data correlation matrix A c . 

25 

39. The method of claim 35, further comprising calculating the 
whitened version r„ w of the linear convolutive environment model r n according 
to: 

r£ =W(H 0 b n + H\b n -\ + n n ) = H 0 b n + H x b n ^ , 
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wherein the linear convolutive environment model r n is represented according 
to: 

r n =n o b n +U 1 b n _ 1 +n n 
where b n and b n _ x are the K-d vectors of current and previous symbol for all 
the K users, H 0 and I^are GxK mixing matrices with the structure 

H o=L H o,o H o,i h o,k1 h i =L H i,o H w ••• H 1K J such that 

H o,k =^ e o2L h iZu , H Jit =V^i2l^ H , and s 0 e 1 represent the energy of the 

1=0 1=0 

current and the previous symbol respectively. 

40. An adaptive detection method, comprising: 

introducing an adaptive weighting matrix into a RAKE structure, 
wherein the matrix is adaptively estimated using at least one of Principal 
Component Analysis (PCA) computational techniques or static Blind Source 
Recovery (BSR) computational techniques based on Independent Component 
Analysis (ICA). 

41. The method of claim 40, further comprising employing a closed 
form RAKE structure according to: 

t \spWH H r forDS-CDMASystems 

i,RAKE-ICAIPCA-\ s H c HwHH ? f orW CDMA Systems 

where w = diag[A A • • • A ] , and A is the matrix. 

42. The method of claim 41, further comprising adapting the matrix 
according to natural gradient update laws. 

43. The method of claim 42, further comprising adaptively estimating 
the matrix using natural gradient update laws according to: 

A(k + 1) = A(k) + ?? k AA(k) 

where 
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AA(k) = 



(I-(p(y(k))y(k) H )A(k) for static BSR (or ICA) 
{ (I-y(k)y(k) H )A(k) for PCA 



and is a nonlinear score function which depends on an underlying 
distribution structure of involved signals. 

44. The method of claim 43, further comprising employing a score 
function according to: 

(Pi (yi ) = V iyi - a t (tanh(# Refe }) + tanh(# Im{;y t })) . 



45. The method of claim 40, wherein a channel estimate H is at least 
one of not available and changes dynamically, the method further comprising 
estimating a detector without using the channel estimate, such that the detector 
structure reduces to Matched Filter BSR/PCA according to: 

% | StfWr for DS-CDMA Systems 

i,MF-BSR/PCA~^ s p c H Wr for WCDMA Systems 
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